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Abstract. 

Spectral properties of the Hamiltonian function which characterizes a trapped ion 
are investigated. In order to study semiclassical dynamics of trapped ions, coherent 
Jlj ' state orbits are introduced as sub-manifolds of the quantum state space, with the 

S . Kahler structure induced by the transition probability. The time dependent variational 

^H' principle is applied on coherent states' orbits. The Hamilton equations of motion on 

Kahler manifolds of the type of classical phase spaces naturally arise. The associated 
classical Hamiltonian is obtained from the expected values on symplectic coherent 
^ ' states of the quantum Hamiltonian. Spectral information is thus coded within the 

^«0 , phase portrait. We deal with the bosonic realization of the Lie algebra of the SU(1,1) 

group, which we particularize for the case of an ion confined in a combined, Paul and 



in 

oo 



f~^ ■ Penning trap. This formalism can be applied to Hamiltonians which are nonlinear 



l> 



in the infinitesimal generators of a dynamical symmetry group, such as the case of 
(3 ' ioiis confined in electrodynamic traps. Discrete quasienergy spectra are obtained 

and the corresponding quasienergy states are explicitly realized as coherent states 
parameterized by the stable solutions of the corresponding classical equations of 
motion. A correspondence between quantum and classical stability domains is thus 
K^ , established, using the Husimi representation. 

^\ 

PACS numbers: 02.20.Sv, 03.65. Ca, 03.65. Ge, 03.65.Fd, 03.65. Sq, 37.10.Gh, 37.10.Ty, 
37.30-Fi 



Time dependent variational principle 2 

1. Introduction. Classical-quantum correspondence 

The development of quantum mechanics and its concepts left it with a heavy legacy of 
classical concepts. The most notable among them is Bohr's correspondence principle, 
which states that in the limit of large quantum numbers there are similarities between 
classical and quantum dynamics. In other words, the results from classical mechanics are 
macro scopically correct and may be regarded as limit cases of the quantum mechanics 
results, when quantum physics discontinuities are negligible. This principle is still used 
as an intuitive guide for finding quantum properties which are similar to known classical 
laws. Some of these analogies are intriguing, considering the fundamental differences 
in the mathematical formalism underlying the two theories: quantum mechanics uses a 
separable Hilbert space with a unitary inner product, while classical mechanics is based 
on a continuous phase space endowed with a symplectic structure [H E]. Nevertheless, 
in certain situations it might represent a weak correspondence between classical and 
quantum concepts. The concept of observables in quantum mechanics and in classical 
probability theory has drawn a lot of interest. It was shown that, by means of injective 
statistical maps, quantum mechanics can to a certain extent be reformulated in classical 
terms [3]. 

Pure and mixed states of a n-dimensional quantum entity can be represented 
as points of a subset of a n^ dimensional real space |1]. Mathematically, a pure 
quantum state is typically represented by a vector in a Hilbert space. Mixed quantum 
states are characterized by means of density matrices. Both mixed and pure state are 
intensively investigated, due to a very large interest which they present for modern 
physics. The expressions coherent and incoherent superposition of quantum states are 
also used in order to help distinguish between pure and mixed states. There exist 
superpositions of quantum states which exhibit many features similar to those of their 
classical counterparts: they are the so-called coherent states, investigated by Schrodinger 
in 1926 13 and then rediscovered by Klauder [6l [7] , Glauber |8l |9] and Sudarshan [10] 
at the beginning of the 1960s. Glauber coined the term 'coherent states' in 1963, in 
the context of quantum optics [H]. Such states are superpositions of Fock states of 
the quantized electromagnetic field that, below the limit of a complex factor, are not 
modified by the action of the photon annihilation operator [12]. Since then, coherent 
states had a huge impact on quantum physics [131 US [HI IISI 113 [HI [IS] with widespread 
applications such as nuclear, atomic and condensed matter physics [20 | |2T | [22 | I231I2H 125] . 
quantum field theory [26], problems of quantization and dequantization [211 ISHl 127] , 
non-commutative quantum mechanics [28] which tries to unify quantum mechanics and 
gravity, path integrals, signal analysis and recently, quantum information processing 
(QIP) [231 [21 [251 |29l |30l [31] using the entanglement feature [E] [231 [261 |32l |33]. 
Moreover, there are a lot of quantum states of interest for modern physics and not only 
[M] . such as Fock states, Schrodinger-cat states [T8l[32] or squeezed states [35 ] 136 ] [37 1 138] . 
A matter of interest is the engineering of quantum states showing negative parts in 
their Wigner functions, which is a signature for non-classical properties of the state 
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that can not be explained using a classical phase-space distribution. All these represent 
problems of ever increasing interest, because they discuss the relation between classical 
and quantum physics, of the blurred separation line between the two of them. 

A charged particle driven by a time-dependent perturbation in a quantum system is 
a non-trivial fundamental issue, with numerous applications such as nonlinear dynamics 
(quantum chaos) [39], quantum simulation of the Dirac equation and implicitly quantum 
field theories |10lllT], cavity QED [l2lll3], quantum non- locality [H], quantum transport 
[I5] . quantum optics [191 SHI 113 HE], superpositions of arbitrary quantum states [19] 
or quantum information processing [32l [331 ISQl [Ml [52l [53]. The paper investigates 
quantum dynamics for systems with dynamical symmetry, such as is the case of ions 
confined in electro dynamic traps or combined (Paul and Penning) traps \5i\. The 
coherent states' orbits are introduced as sub-manifolds of the quantum state space, 
with a Kahler structure induced by the transition probability. An algorithm was found, 
which associates a classical Hamiltonian to the quantum Hamiltonian describing the ion. 
The classical Hamiltonian implicitly contains the spectral information of the quantum 
system. The coherent states formalism introduced in [551 ISS] and developed in [27] is 
used, and the time dependent variational principle (TDVP) [571 EHl [59], in order to 
characterize the semiclassical dynamics. 

The results have been applied to the case of an ion confined in a quadrupole 
electromagnetic trap with anharmonicities, in order to study semiclassical dynamics. 
The coherent states formalism for dynamical groups and the time dependent variational 
principle (TDVP) are applied with an aim to investigate semiclassical behaviour of 
trapped ions. The explicit expression of the quantum Hamiltonian associated to a 
charged particle is found, as a function of the Lie algebra generators for the radial 
and axial symplectic groups. The Schrodinger equation solution also results. The 
Hamiltonian function is then particularized to the case of combined, quadrupole and 
octopole traps. 

The paper is organized as follows: Section [2] reviews some of the definitions from 
the geometry of coherent states. In Section [3] the time dependent variational principle 
is applied to a manifold of test vectors in paragraph 13. 1[ The Lagrange function 
is introduced and the Hamilton equations of motion result for the system we have 
investigated in Section 13.21 In Section 13.31 we introduce the coherent states built over 
the weight states of the discrete positive series of the symplectic group Sp (2, M) and 
the Lie algebra generators. A model is proposed which enables the explicit calculus of 
the Husimi function for every algebraic model with a dynamical group. The equations 
of motion in phase space are inferred. The model we suggest is applied in Section [4] 
to a combined, Paul and Penning trap, with octupole anharmonicity. The Lie algebra 
generators for the radial and axial symplectic group Sp (2, R)^ ^ are introduced in Section 
14.11 together with the Hamiltonian of the system. The symplectic coherent states 
representation is introduced in Section 14.31 The classical Hamiltonian for the system 
we investigated is finally obtained. Section [5] is dedicated to a discussion of the results. 



Time dependent variational principle 4 

2. Coherent states for dynamical symmetry groups 

2.1. Notions of differential geometry 

Euclidian space is a basic symplectic manifold, not only because it is simple to 
characterize, but mainly to the fact that it provides a local model for every symplectic 
manifold. The Euclidean space can be treated as a semi-local symplectic manifold. 
Hence we can submit global nonlinear problems on it. The phase space is a space 
endowed with a tensor, which gives it the geometric structure. It is called a symplectic 
form. Symplectic manifolds are even-dimensional. They are studied using symplectic 
geometry or symplectic topology [601 El], which represents a branch of differential 
geometry (topology) [62] . Symplectic manifolds arise naturally in abstract formulations 
of classical mechanics and analytical mechanics as the cotangent bundles of manifolds, 
e.g. in the Hamiltonian formulation of classical mechanics [601 ES]- The subject of 
symplectic topology is the global structure of a symplectic manifold and the behaviour 
of symplectomorphisms which are far from identity. These aspects are of interest both 
in Euclidean space with its associated standard linear symplectic structure and on 
general symplectic manifolds [611 E2]- A Kahler manifold is a symplectic manifold 
with an integrable almost complex structure. A consequence of the integrability 
theorem is the fact that every Riemann surface carries a Kahler structure [6^. Most 
symplectic manifolds are not Kahler, which means they do not possess an integrable 
complex structure compatible with the symplectic form. Gromov [651 ES] made the 
important observation that symplectic manifolds admit an abundance of compatible 
almost complex structures, so that they satisfy all the axioms for a Kahler manifold 
except the requirement that the transition functions be holomorphic. Any real valued 
differential function if on a symplectic manifold can serve as an energy function or 
Hamiltonian. 

2.2. Coherent states geometry 

Further on, we will recall some definitions and constructions from the geometry of 
coherent states [131 El EH], while introducing the elementary quantum systems (in 
the Wigner sense), which represent the most extended class of models describing both 
analitically and numerically the energetic spectra and the solutions of the Schrodinger 
equation. 

Let G be a Lie group acting on a manifold M (respectively, a complex Lie group 
acting on a complex manifold M). Then for every point m G M we define its orbit by 
Om = Gm = {9-i^\g € G} and stabilizer by [67] 

G^ = {ge G\g.m = m} (1) 

We will denote by H the separable complex Hilbert space of the state vectors 
(quantum states). The state space in quantum mechanics will be denoted by V{T-L). If 
V' G "H is a state vector, then ifj = {XiplX G C} represents the state associated to ip- 
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v{n) = \^ ^pen\{o}\ , (2) 



is the state space, which is also a projective complex space, ip is an unidimensional, 
complex linear space. If using normalized states {\\ip\\ = 1), we can introduce the states 
as radii and we have -^ = < e**^ ■ V^ y? G M L In the following, we will consider U as an 
irreducible unitary representation of the Lie group Q in "H. This means that V({g) is a 
unitary operator on Ji for \/g G Q, such as 

U{g-g')=U{g)-U{g'), ^g,g'eg 

U{g-') = [U{g)]-'=U+{g) . 

where U^ is the adjoint of U. This representation induces an action of the group Q on 
ViTi), defined as gip = U {g)ip for Wg E Q and tp G "HyjO}. Then Qq is the stationary 
group of the state ip, defined as 



Go = [g g-'^ = i' , g ^Gj 



The orbit G ■ i^ = \g ■ i' fi'^^fis called a coherent state orbit, while any state g ■ ip 

is called a coherent state for the group Q. We will denote Ox = Q ■ ip. The coherent 
vectors which determine the coherent states, can be parametrized with the coordinates 
of the quotient space A^ = Q /Qo- For a large class of symmetry groups Q, M. presents 
a structure of classical phase space (symplectic manifold). Thus it is convenient to 
parametrize the try functions with points from M.. There exists a bijection between the 
points from M., considered as classical states, and the coherent states. 

The triplet {1-L,U,Q) is defined as an elementary quantum system in the Wigner 
sense or a quantum system with symmetry, where "H is a separable Hilbert space, Q is 
a group and U stands for an irreducible unitary representation. Coherent states can be 
associated with such systems. We will consider that the orbit Oj, is a sub-manifold of 
the space V {1-1) with a structure of Kahler space induced by the transition probability. 
Coherent state Kahler orbits are obtined for Heisenberg groups, compact groups and 
semisimple groups, with representations belonging to the holomorphic discrete series. 
The coherent states considered here have important applications in quantum physics and 
quantum optics, especially when referring to the physics of trapped atomic particles. 

Thus a family of coherent states provides the opportunity of investigating the 
quantum world using the methods of study applied to the classical world. This feature 
allowed Glauber and others to treat a quantized boson or fermion field as a classical 
field, with an aim to compute correlation functions and other quantities relevant for 
statistical physics, such as partition functions. Hence, the system can be characterized 
by means of the associated classical trajectories. 
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2.3. SU(1,1) group in quantum optics 

We now turn our attention to elements of the group theory, especially those related to 
the SU(1,1) group, which presents a great interest for this paper. A. O. Barut and L. 
Girardello [68], then A. M. Perelomov [131 IH] have brought in the SU(1,1) group to 
quantum optics. They have constructed the coherent states of the SU(1,1) group in 
different space representations: Hilbert space and unit disk space or Lobachevian space 
[111169]. We will focus on geometrical aspects of coherent states and generalized coherent 
ones, based on Lie algebras, especially on the su(l, 1) Lie algebra [TOl [TH [721 [73]. The 
topological space that parametrizes the elements of a Lie group is a manifold. We can 
ascertain that a Lie group is a finite dimensional manifold. 

Generally, group methods are used in order to find degenerate states of a given 
state in quantum mechanics [711 EH] . The Lie algebra of the generalized special unitary 
group SU(1,1) is widely used in quantum optics. The SU(1,1) coherent states arise in 
the case of generation of boson pairs with zero spins in homogeneous alternating electric 
field or in the gravitation field of the expanding Universe [76] . In particular Wodkiewicz 
and Eberly [77] have discussed the role of the generalized coherent states of Perelomov 
[131 [H] associated with the Lie algebra of the SU(1,1) group in connection with variance 
reduction. Applications of this group are tightly connected with squeezed states of hght. 
Coherence-preserving Hamiltonians associated with the SU(1,1) generalized coherent 
states have been studied by Gerry [791 [HO]- It should be emphasized that the SU(1,1) 
generalized coherent states in [13] and [77] are special cases of the two-photon coherent 
states discussed by Yuen [3H]. This fact has been used by Gerry [THl [HD] who applied 
the SU(1,1) formulation of the two-photon coherent states to the issue of the interaction 
of squeezed light with a nonlinear non-absorbing medium modelled as an anharmonic 
oscillator. The coherent states have been especially useful in quantum optics and 
quantum information processing (QIP) [El [23l [24]. Each mode of the electromagnetic 
field may be formally described as a harmonic oscillator, and different quantum states 
of the oscillator correspond to different states of the field. The field from a single- 
mode laser operating far enough above threshold can be described for many purposes 
as a coherent state; it differs from a coherent state in that its phase drifts randomly. 
Coherent states also play an important role in mathematical physics [121 El EHl [73] . In 
the present paper we will study the bosonic realization of the Lie algebra for the SU(1,1) 
group. 

The 5^7(1, 1) group [12] consists of unimodular matrices of 2 x 2 dimension with 
unit determinant. It is also isomorphic to the following noncompact groups |81j : 



50(2, 1) ^ 5f/(l, 1) ^ SL{2, M) ^ Sp{2, M). (3) 

For example, we can construct the three-dimensional noncompact group S'f/(1, 1) 
by taking the intersection of the unitary group f/(l, 1) with SL{2; C). All usual groups 
of linear algebra such as G'L(n,M), 5L(n,M), 0(n,M), U{n), SO{n,M), SU{n), Sp{n,R) 
are real or complex Lie groups [67j. 
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3. The Variational Principle 

3.1. The Variational Principle. Test states 

The equations of motion in classical mechanics arise as solutions of variational problems. 
An example would be Fermat's principle of least time, stating that light propagates 
between two points by a path which takes the shortest possible amount of time. In a 
manner similar to light, all systems possessing only kinetic energy move along geodesies, 
which are paths minimizing energy. A general mechanical system possesses both kinetic 
and potential energy. In such case the quantity which is minimized is the mean value of 
the difference between the kinetic and potential energy. This less intuitive quantity is 
called the action [61] . The variational principle results from the definition of an action 

integral S = J L {ip) dt , where the real Lagrange function of this variational problem is 

defined as ^ El 



LW ^ 



\H\ih) — Qm ' 



^ 



(4) 



{i}\i))Y ' ' ' \dt 

where KH is the quantum Hamiltonian of the system. In eq. (jl]) ■?/' is a vector in the 
Hilbert space H, for any moment of time t, and it belongs to the domain of the self- 
adjoint quasienergy operator Kit) = H — id/dt. The variational principle for S applied 
on the entire Hilbert space 1-i, yields to the Schrodinger equation ^7\ 15^. 

By minimizing the action {5S = 0), it can be shown that the Schrodinger equation 
is rigurously obtained from the variational principle. The natural symplectic structure 
induced by the transition probability between states for unitary transformations, 
achieved by the identification of the symplectic form with the imaginary part of the 
scalar product, enables approaching quantum mechanics issues through the formalism of 
the classical mechanics on (infinitely dimensional) symplectic manifolds (phase spaces), 
in particular Kahler manifolds. 

We apply the variational principle to a manifold M. of test vectors, parametrized by 
the points of a finite 2n dimensional phase space M., with a Kahler manifold structure. 
In agreement with the Darboux theorem, the symplectic manifold M. admits canonical 
local coordinates and associated complex local coordinates. Generally the complex 
structure is not global, but local complex parametrization is preferred. In case of 
elementary quantum systems with dynamical symmetry groups which admit coherent 
states, the complex structure is global and M. represents a Kahler manifold. 

3.2. The Lagrange function. Hamilton equations of motion 

We will consider z = {zi,Z2, ■ ■ ■ ,Zn) G (9 as a system of complex canonical local 
coordinates in Ai, where O is an open set from C". A^ is a manifold of dimension 
2n. The (global) symplectic structure and the (local) Kahler structure are induced on 
the variety of test vectors Ai. We choose a family of vectors ip {^) G A^ with z & O, 
holomorphic in z, such as dip {z)/dz* = 0, i = 1, . . . ,n . The Lagrange function in the 
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complex parametrization can be expressed as 



L (z, z 



{i:{z*)\H\^{z)) 



E^A^(^*) 



dip (z) 



dzi 



— z- 



dip (z*) 



dz-i 



^W) - 



{ij{z*)\ij{z)) • 
It is convenient to introduce the following notations [57] : 

N{z,z) = {'iP{z)\,p{z)) , H,i{z.z) = ^^^-^-^^y-^-^-^ , (6) 

where H is the quantum Hamiltonian of the system, while H^ stands for the classical 
Hamiltonian. Then the Lagrange function described by eq. ([5]) becomes 



L{z, z*) = -{zV,- z*V.,)\YiN {z,z*)- H,i{z,z*) , 

n 

where z -V z = ^ Zidz^ ■ Hence the variation of the action is 



(7) 






(8) 



6S = / < z > - — ;r— z-iSzt - i >^ z*6zk - SHci 

j,k=i ■' "- j",; 

From the condition of extreme (55* = 0, using independent variables, we obtain a system 
of equations such as : 



v^ 92 In AT 



dz*dzk^^ 






dz 



jr[ dzjdzl 



dzl 



We now introduce the matrix of the symplectic structure on M. 

92 In N 



n 



\.^jk)l<j^k<n 



^jk 



dzjdzl 



and Eqs. (|9]) change accordingly to 



E 



UJkjZ 



■<)^3 



dzk 



i 2^ (^jkZj 



dH, 



d 



dzl 



(9) 



(10) 



(11) 



where u% = Ukj- Hence we can ascertain that the matrix Q = Uij is Hermitian. The 
Poisson brackets for the f,g^ C°° (A^) functions, smooth on Ai, can be introduced as: 



n , 

{f,9} = i Yl (^i^ 



df dg dg df 

— A, 



dzjdzl "^^dzjdzl 



dz I dz* \dz j dz 



(12) 



,-1 



where \jk = —iiVt*) ^ . Thus the variational principle for S applied on M. leads to the 
following equations of motion, which are classical [57]: 



f ^{.,i..},f M.%i..} 



(13) 



(5) 
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where 

*=''-'^- (■v.MiV'W) ■ <"* 

This paper deals with the study of classical and semiclassical nonlinear dynamical 
systems. It focuses on similarities and differences found when comparing the dynamical 
behaviour of Hamiltonian systems described by classical theory in respect with those 
described by the quantum theory |82l [83] . An algorithm has resulted, through which 
an energy function is associated to the quantum Hamiltonian. This function represents 
a classical type Hamiltonian, whose values are precisely the expected values of the 
quantum Hamiltonian |12j. The energy function determines nonlinear equations of 
motion which can be approached using the theory of differential dynamical systems. 
It is suggested that this algorithm could be used in order to establish a correspondence 
between semiclassical and classical dynamics. 

The phase space A^ is a symplectic manifold described by eq. ( fTOj) . with Poisson 
brackets determined by eq. ( TT3l) . Eqs. ([9]) are the Hamilton equations which describe the 
Hamiltonian function Hd on Ai. Thus Hd {z, z*) is the expected value of the quantum 
Hamiltonian in the state represented hj iIj{z) G Ai. Consequently Hd is considered the 
classical Hamiltonian associated to the quantum Hamiltonian H. This association is 
usually called dequantization [HI] and it represents a set of rules which turn quantum 
mechanics into classical. The solutions of the equations of motion for the classical 
Hamiltonian Hd contain the same amount of spectral information as the solutions of 
the Schrodinger equation for the quantum Hamiltonian iJ, for many of the dynamical 
systems of interest. One of the conditions to be fulfilled is to choose A^ as a manifold 
of coherent states. 

The group of linear canonical transformations of a dynamical system with n degrees 
of freedom is the symplectic group Sp {2n, M). An electrically charged particle confined 
within an electromagnetic trap can be explicitly described using coherent states, for a 
subgroup Q of the symplectic group 5p (6, M) [271 129]. Symplectic transformations in 
Sp (2n, M) leave invariant the canonical form of the classical Hamiltonian equations of 
motion [611 EZ] • 

3.3. Symplectic coherent states. Algebraic models. Equations of motion 

Further on we will introduce the coherent states built over the weight states of the 
discrete positive series of the symplectic group Sp (2, M), which is also a real Lie group 
[271 l55l l56l [85] . We will consider a unitary representation U of the group Sp (2, M) in 
the Hilbert space "H. We will denote by Q the group of unitary operators U (g) with 
g E Sp (2, M) and by g the Lie algebra for this group (in fact, the SU{1, 1) Lie algebra). 
The algebra is determined by three generators Ko,Ki, and K2. Further on we may 
define the ladder operators K± = Ki ± iK2 so that the Lie algebra may be cast into 
canonical form [3^ 



[Ko,K^] = ±K^, [K_,K+] = 2Ko. (15) 
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The Lie algebra sp(2,M) consists of all the real linear combinations of the operators 
iKQ,iKi, and iK2. We will denote by U^ a unitary irreducible subrepresentation of 
the U representation in the Hilbert subspace Hk of the "H space, where k stands for 
the Bargmann index [86] . The discrete series unitary representation of the Lie algebra 
considered are labelled by the eigenstates of the Casimir operator C [56] 

C2 = Kl - Kl - Kl , (16) 

with eigenvalues k{k — l) [Til EZl [88]. We will consider only irreducible unitary 
representations of the group Sp{2,'R), corresponding to the discrete positive series 

v+{k) [MEniEl]. 

Coherent and squeezed states in Quantum Optics are concerned with the Positive 
Discrete Series D^{k). The canonical ortonormal basis of the Tik space consists of 
the \m,k); m = 0,1,..., vectors. For the representations of interest, the states 
\m, k) diagonalize the compact operator Kq. The operators Kj^ and K_ are Hermitian 
conjugates of each other and act as raising and lowering operators of the quantum 
number m [311 ESI E] : 

K_^ |m, k) = [{m + 1) (m + 2k)]^^'^ \m + l,k) , 

K_ \m, k) = [m{m + 2k- l)]^'"^ \m -l,k) , (17) 

Kq \m, k) = {m + k) \m, k) . 

The coherent states for the symplectic group 5p(2,M) are defined as |34] : 

\z^m^k) = U {z)\m,k) , (18) 

Using the disentagling theorem for the S'f/(1, 1) Lie algebra [771 [78] the unitary operator 
U{z) can be expressed as 

U {z) = exp {zK+) exp {l3Ko) exp (-zK^) , (19) 

where \z\ < 1 and f3 = In {1 — zz). The |z, m, fc) states are called symplectic coherent 
states. The phase space M. whose points parametrize the coherent states is the unitary 
disk \z\ < 1, endowed with the Lobacevski metrics ds'^ = 4 (1 — zz)~ dzdz. For m = 
we infer the geometrical construction of Perelomov [H]. Eq. ( ITSll remains valid for 
coherent states with m > 0, realized by applying the unitary operators U {z) E Q on 
the nondominant weighting vectors \k,m) ] m = 1,2, . . . . 

Using the Baker- Campbell-Hausdorff formula [671 [HH] and eq. ( !T8|) . we infer 
[Ml [561 [90] 

U-' (z) KoU (z) = (1 - zz)-^ [zK^ + (1 + zz) Kq + zK+] , 

U-'^{z)K+U{z) = {l-zzy'^[K+ + 2zKo + z^K_] , (20) 

U-^ (z) K^U {z) = (1 - zzY^ [K^ + 2zKq + z^K+] . 

We will denote by fi^ = Kq + eKi, where e = ±1. Then 

i.-.(,)n,t'(,)=(i±££l(l±££)£,, (21) 

Ee = \{l- zz)-^ [2 (1 + 2) (1 + 2) Ko + £ (1 + zf K_ + £ (1 + zf A'+] (22) 
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From eqs. ([20]) - <^, we infer 



(z, k, m\Q'^\z, k,m) 



[l + ez)(l+ez) 



{k,m\E^\k,m) . (23) 



1 — zz 

Eq. f l23|) enables the explicit calculus of the Husimi function for every algebraic model 
with a dynamical group Q. If H is the Hamiltonian for such a model, expressed as a 
polynomial in the operators Kq, Ki, and K2, we observe that the Hamiltonian function 
Hd (z) = {z, k, m\H\z, k, m) leads to the following equation of motion in the phase space 
M : 

z = {z,Hd} . (24) 

The brackets { , } stand for the Poisson brackets associated with the symplectic 
manifold M., defined as [5 



/f 1 _ (1-^^)^ (dfdg dfdg\ 

^•^'^^ 2i{k + m)\dzdz dzdzj ' ^ ' 

where the phase space classical observables / and g are smooth and real functions of 
^ez and '^mz. When m = 0, from the TDVP applied to the coherent states \z, 0, k) 
for the Hamiltonian H we infer Eqs. ( 12^ . which can also be expressed as 

_ (1 - z-zf dHd ^26) 



2ik dz 
We now introduce the complex variables ^ and rj defined as 

i = (^ + ^)(^ + ^) , = (1-^)(1-^) . (27) 

1 — zz ' 1 — zz 

Then eqs. (I26l) become 



^ _ 2 z- z dHci . _ 2 z- z dHci , 

ik 1 — zz drj ' ikl — zz dC, 

which are exactly the equations of motion in the phase space for the system we 
considered. The group of the linear canonical transformations of a dynamical system 
system with n degrees of freedom, is the symplectic group Sp(2n,M.) [HI]. In case of 
dynamical symmetry groups, the Schrodinger equation solutions associated to linear 
Hamiltonians are given by coherent vectors multiplied by geometrical (Berry) phase 
factors [SS]. The associated classical Hamiltonians result from the expected values of 
the quantum Hamiltonians on coherent symplectic states. Thus spectral information 
is coded into the phase portrait. The formalism presented here allows to explicitly 
construct bases and systems of symplectic coherent states for the study of trapped ion 
systems. 

4. Semiclassical dynamics of an ion in combined traps with axial symmetry 

4.I. Algebraic models for nonlinear axial traps 

We have reviewed the properties of the coherent states and their exceptional importance 
for physics and especially for quantum optics. Coherent states for charged particles 
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confined in electrodynamic traps can be built using the group theory [211 [S2]- We 
will further investigate the semiclassical dynamics of an ion trapped in a nonlinear 
electrodynamic trap with cylindrical symmetry, using the symplectic state approach 
and coherent state formalism developed in [HI |271 [291 ES] ) presented in Sections 13.21 
and I3.3[ The coherent state formalism is applied in order to study the semiclassical 
behaviour of trapped ions [271 ISll [551 [56]. We will introduce the Lie algebra generators 
of the axial symplectic group Sp (2, M)^ [271 ISS] 



Ko. 



la 






z^± 



Pz 



M'^LO?. 



Ko„, = -^ 



i 


\ d 1 




2z ^ +1 


Ah 


I 9z \ 



and the Lie algebra generators of the radial symplectic group Sp (2, 
eigenvalue of the orbital angular momentum hi 



K, 



0,lr 



2hUJr 



2M 



p;± 



Mul 



-P' + 



2M 



Kor = -r- 



Ah 



2p ^ + 1 

op 



(29) 
for a fixed 

(30) 



We have denoted by uJa/2T: the axial frequency and by ujrl2T: the radial frequency 



respectively, of the quadrupole trap potential. The Koj,Kij and 7^2^ 



i[Kij,Koj],j 



a, r operators satisfy the commutation relations for the Lie algebra of the symplectic 
group Sp (2, R) 

[Koj,Kij] = iK2j, [K2j,Koj] = iKij, [K2j,Kij] = iKoj. (31) 

Then the quantum Hamiltonian of a particle with electrical charge Q, mass M 
and orbital angular momentum hi, located in the quadrupole trap with octopole 
anharmonicity, can be expressed as [90] 



Moo a iOr 



U 



2Q 



4 + M "^ ^^^ ' ^^°^' + ^'"'^ ~ 



— hl + hhJa [Koa 



Kia) + hWr {Kor - i^lr) + C2A (t) H^ (p, z) 



(32) 



where Uc stands for the cyclotronic frequency and H4 is a harmonic polynomial of degree 
2 homogeneous in the variables z^ and p^ = x^ + y'^. These variables can be expressed 
as a function of the generators 



Ah 

MUr 



(Koa + K,a) 



P 



Ah 

MUr 



{Kor + K 



Ir) 



(33) 



The coupling constant C2 and the time dependent function A (t) are characteristic to 
the type of trap and will be explained later. 

Eq. fl32|) results from the general expression of the Hamiltonian associated to a 
particle which undergoes the influence of an electromagnetic field, characterized by an 
electrical potential $ and a magnetic induction B [191 Ell EHl [85] : 

2 



H 



2M 



1 



p 



QB xf 



g$ . 



(34) 
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The inomentum operator for the particle is p = —ihV, while the magnetic field is 
considered to be axial: B = (0, 0, Bq) and Uc = QBq\M. We can introduce the (axial) 
angular momentum operator on the z axis, defined as Lz = xpy — ypx- Then Eq. ( IMj) 
becomes 

ff = ^W.+P?+p;) + sij«=i5o^r-|§i. + 0*. (35) 

Because the axial angular momentum operator L^ commutes both with the symplectic 
groups generators as well as with the Hamiltonian H, the study of the quantum system 
can be restricted to a subspace of the Hilbert space for which L^ = hll, where / is the 
orbital quantum number. The Hamiltonian reduced to this subspace can be expressed 
as 

ff, = ^(p^+p; + p3 + 5i^OXP=-||«H-0*. (36) 

Choosing a system of cylindrical coordinates p and 6, the expressions for the radial 
coordinates are x = pcos6 and y = psinO. It follows that 

2,2 ^,ja^ 1 d-' , id\ , ^^ a 



We choose \Ef = e*'^x as a solution of the Schrodinger equation, where x depends only 
on p and z. Then L^^ = fd^ and x is a solution of the Schrodinger equation for the 
Hamiltonian 

which results from Eqs. ( l36l) and ( 1371) . 

We will now choose a potential with axial symmetry, such as 

^{f,t) = A{t)g{p,z) , (39) 

where A if) is a time-periodical function of period T = 2tx/VL, while g is a. function of p^ 
and z^ which satisfies 

9 (P, ^) = X^ CkH2k (p, z) . (40) 

k>l 

The H2k polynomials are harmonic polynomials of degree k in p^ and z"^. The potential 
described by Eq. ( l39l) is characteristic for a trap with axial symmetry and with symmetry 
in respect to the radial plane xOy. For a Penning trap A is constant in time, but for 
most of the radiofrequency traps we choose A{t) = Uq + Vq cos Vtt. In case of a harmonic 
potential (c^ = for /c > 1) we have 

1 1 

^(P'^) = ;:^-^(p'-2^')' c2 = --j-— 2, (41) 

' "T ^^0 ' "T ^^0 

where tq and zq are the semiaxes of the combined quadrupole trap we have considered. 
We remind that a combined trap consists of a Penning trap and a superimposed Paul 
trap. The particular case of the ideal Paul trap is obtained for i? = 0. The quadrupole 
potential is produced by hyperbolic electrodes. Real traps are different from ideal ones. 



Time dependent variational principle 14 

Due to mechanical imperfections (holes drilled in the endcap electrodes or in the ring 
electrode, the latter causing a radial asymmetry of the potential), the trap potential is 
not a pure quadrupole potential. The potential is then expanded in a power series and 
higher order terms have to be considered. Due to the trap symmetry we have chosen, the 
odd terms of higher order from the series expansion of the potential vanish and we are 
left only with the even terms. In case of a quadrupole trap with octopole anharmonicity, 
the electric potential can be expressed as 



g (p, z) = C1H2 (p, z) + C2Hi (p, z) , 




(42) 


H2{p.z)=2z^-p\ H,{p,z) = 8z^- 


-24^V + 3p^ 


(43) 



where H2 is the pure quadrupole term and H^ is the octopole term. Under such case, 
the Hamiltonian Hi is given by Eq. (!32l) . The next anharmonical contribution is given 
by the polynomial 

(44) 



He = 


: 16Z^ - 


120zV + 90z^p^ - 


-5/. 






By denoting 














Kr = 


^f 2Qc2A{t) 


, K,-- 


= 4Qc2A{t) , 


Hamiltonian given by Eq. (j38|) can be expressed 


as 




Hi = 


2M 


ydp'^ p2 pdp) 


2Mdz^ 


-f« + 


-t 


2 , Kc 

^+ 2 


'-z^ + QA{t)P{p\ 


-^)> 







(45) 
where the anharmonical part is defined as 

P{p\z^)=Y,CkH2k{p,z) . (46) 

fc>2 

In the following we will study the Schrodinger equation 

dx 

where Hi is given by Eq. ( HSl) . with particularizations for combined quadrupole and 
octopole traps. In agreement with Eqs. ( 132|) . (133|) and ( H5|) . The Hamiltonian Hi 
describes an algebraic model when P{p'^,z'^) is a polynomial function. This model is 
linear for quadrupole traps [P is a linear combination between p^ and 2;^) [90] . 

4-2. Dynamical symmetries for an ion in an axial trap 

The algebraic model described by the Hamiltonian Hi can be characterized by the 
dynamical symmetry group Q = Qa ® Qn realized as a direct product between the 
axial symplectic group Qa = Sp (2, M)^ and the radial symplectic group Qr = Sp (2, M)^ 
ESI [93]. 
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From Eqs. ( 130 p we obtain the commutation relations 

for the Lie algebra Qa of the group Qa- The Casimir operator of the Lie algebra Qa 
determines the Bargmann indexes /c, and it can be expressed as [551 EHl [71] 

C,^ = Kl- Kl- Kl = -^I = k{k-l)I , (48 ) 

where / is the unitary operator. From Eq. fllSl) we infer the values of the two Bargmann 
indexes which characterize the axial motion 

K+ = \ , K^ = \- (49) 

From Eqs. (!30|) we obtain the commutation relations 

[KQr.Ki,] = iK2r , [K2r,Kor] = iKir , [-R'lr, -?^2r] = ^-^Or , (50) 

for the Lie algebra Qr of the group Q^- The Casimir operator of the Lie algebra q^ is 
[32 [561 [90] 

C2. = Kl - Kl - Kl = ^-^ I = k{k-l)I . (51) 

From Eq. (IHTj) we infer the Bargmann index /:;,, = (/ + 1) /2, where / represents the 
quantum orbital number. Eq. fj^5|) can be now expressed as 



Hi = hUr {Kor - Kir) + 2/lWa (-^^Oa " -ft'la) " y^/ 

Consequently, the study of the Hamiltonian for an ion confined in a quadrupole trap 
(P = 0) has been restricted to the study of a linear model for the dynamic group Q, for 
which the irreducible unitary representations are characterized by the Bargmann indexes 
ka and kj., fixed |55]. The solutions of the Schrodinger equation for this elementary 
quantum system will result explicitly by using the coherent state formalism for the 
group Q. This formalism will be also used for algebraic models which describe nonlinear 
traps (with P 7^ 0). 

It can be shown that symplectic coherent states are realized in the trapped ion 
dynamics and a classical Hamiltonian (the Husimi or the energy function) can be 
associated to the quantum Hamiltonian. The Hamiltonian function leads to equations 
of motion in the classical Lobacevski phase space. 

4-3. Symplectic coherent states representation 

According to paragraph 14.21 the Bargmann indexes for the axial motion are given by 
k_ = ^ and k^ = |, while in case of the radial motion the expression is ki = {1 + 1) /2, 
where / stands for the orbital quantum number. The canonical orthonormal basis of the 
Hilbert space for an irreducible unitary representation of the group Qj, of Bargmann 
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index kj, with j = a, r, consists of the vectors (pjm , m = 0,1, . . . and it is characterized 
by the equations 

Koj(f)jm = {kj + m) (j)jm, K^j(j)jj = 0. (53) 

We can introduce the coherent symplectic vectors ipjm (zj), defined as 

^jm (Zj) = Uj {z) (pjm, (54) 

Uj (zj) = exp (zjK+j) exp {PjKqj) exp (—ZjK_j) , (55) 

where Uj (zj) are unitary operators of the representation, /3j = In (1 — ZjZj), and Zj is a 
complex variable from the unit disk \zj\ < 1. The phase space Ai associated to an ion 
confined within an axial symmetry trap represents the product of the axial and radial 
unit disks, realized as a bound domain in C^, made out of all points {za, z^) with l^;^! < 1 
and \zr\ < 1. The coherent states for the dynamical group Q, parametrized with points 
from M., can be expressed as 

In case of quadrupole traps, the symplectic coherent states ^mamr {za, z^) multiplied 
by Berry phase factors are solutions of the Schrodinger equation [55] . 

The energy function associated to the Hamiltonian (l52l) represents a classical 
Hamiltonian Hd on Ai, whose values are exactly the expected values of the Hamiltonian 
Hi on coherent states ^mamr {za, Zr), such as 

Hd = hwr {K + m,) ^^ "/'-^ ^^~ ''^ + 2nwa {ka + m,) ^^ " '"^ ^^ ~ ^"^ 



1_ I ^ y-O. I •■-LLj -| _ 



MuJr 1 — ZrZr MUa 1 — ZaZa 

- -^hl + Hanarm , (57) 



cu 



Hanarm = QA{t) {P {p\z^)) = QA{t)Y,Cu{H2k{p,z)) . (58) 

k>l 

We have denoted by (X) the expected value of the operator X in the coherent state 
^manir {za, Zr) ■ In particular we obtain 

{H^) = 8S2a - 24SirS2r + 3^2r , (59) 

{H^) = 1653a - 120S2aSlr + 90S,aS2r " 5^3. . (60) 

where 

Ojr = Oj [Zr, Kr, TTlr) , Oja = Oj [Za, fca; f^a) i 1 ^ J ^ 3 , 

with 

-:i + z)(i + z) 



Sj {z, k, m) 



1 — zz 



3 



Qj{k,m) , (61) 
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and 

Qi (A;, m) = 2 (A; + m) (62) 

Q2{k,m) =2k{2k + l) + 12km + 6m^ (63) 

Qs (fc, m) =4:k{k + 1) {2k + l)+4mA; (5 + 12A;)+4m^ (15A; + l)+20m^(64) 

The values of the indexes are : 

Tfla, rfly U,i,... , /C(j ~A^ ~A ^ O ' 

We will introduce the complex variables ^j and r^j, j = a,r, defined as 



:i + zj) (1 + zj) 



Zj) [L - Zj, 



Then the classical Hamiltonian (the Husimi function) for an octupole trap in case of the 
pseudopotential approximation can be expressed as 

+ {D,oC + D^iCia + D,^Ul + D,,ea) - Y^l (66) 

If P is a polynomial of p^ and z"^, then Hd is a polynomial of ^^ and C,r- In order to find 
the minimum points of the Husimi function we have to solve the equations 
dHd dHd dHd dHci 

O^a O^r Or]a Orfr 

It can be observed that the Hamiltonian function Hci= {z,k,m\H\z,k,m) 
determines the following equation of motion in the classical Lobacevski phase space 

z = {z, Hd} , (68) 

where { , } stands for the generalized Poisson bracket, defined in Eq. ( !25l) [57j. After 
applying the TDVP to the coherent states \z, 0, k) for the Hamiltonian H we infer Eqs. 

which represents a classical type equation of motion. If P given by Eq. ( H6|l is a 
polynomial in p^ and z"^, then Hd is a polynomial in ^a and C,r- 

The results obtained here can be extended to a system of N particles (ions), 
confined in a quadrupole electromagnetic trap with cylindrical symmetry [561 123 |29] . 
As shown in these references, the Hamiltonian of the center of mass (CM) is similar to 
the Hamiltonian for one particle H, given by Eq. (IMj) . Then the quasienergy states 
and the coherent states are exactly those for the H Hamiltonian, where we have to 
substitute M the mass of one particle with NM. In case when the interaction potential is 
translation invariant and homogeneous of degree —2, we obtain soluble models (e.g., for 
Calogero type potentials). These models enable to explicitly obtain bases and systems 
of symplectic coherent states, in order to investigate ordered systems made of trapped 
ions. 
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5. Conclusions 

The widest class of quantum models which describe both analitically and numerically 
the energy spectra and the solutions of the Schrodinger equation, is made of elementary 
quantum systems (in the Wigner sense). An elementary quantum system is described 
by an irreducible unitary representation W of a dynamical group ^ on a Hilbert space 
Ti. The space of states is the quantum phase space 'P(H), with the symplectic structure 
induced by the transition probability. It is assumed that the representation U admits 
coherent states. Then, coherent states manifolds are orbits of the action of the dynamical 
group on the space of states. The quantum Hamiltonian is a polynomial in the generators 
of the Lie algebra q of the group Q. This property is specific to dynamical groups 
considered as spectrum generating groups. The energy (Husimi) function associated to 
the quantum Hamiltonian if is a classical type Hamiltonian, whose values are precisely 
the expected values of the quantum Hamiltonian on coherent states. 

Thus the quantum dynamical system is considered as a Hamiltonian system over 
the phase space of quantum states, while the corresponding classical system is described 
by the energy (Husimi) function over the phase space Ai, whose points parametrize the 
coherent states. The symplectic structure of the classical phase space and the Hamilton 
equations of motion are explicitly obtained using the variational priciple applied on 
coherent states and from the theory of group representation. 

This paper deals with the study of classical and semiclassical nonlinear dynamical 
systems. It focuses on similarities and differences found when comparing the dynamical 
behaviour of Hamiltonian systems described by classical theory in respect with those 
described by the quantum theory [82]. An algorithm has resulted, through which an 
energy function is associated to the quantum Hamiltonian. This function represents 
a classical type Hamiltonian, whose values are precisely the expected values of the 
quantum Hamiltonian. The energy function determines nonlinear equations of motion 
which can be approached using the theory of differential dynamical systems. It is 
suggested that this algorithm can be used to establish a correspondence between 
semiclassical and classical dynamics. 

Semiclassical dynamics of trapped ions has been investigated by introducing the 
coherent state orbits as sub-manifolds of the quantum states' space with a Kahler 
structure induced by the transition probability. In case of dynamical symmetry groups, 
the solutions of the Schrodinger equation associated to linear Hamiltonians (in the 
infinitesimal generators of a dynamical symmetry group) are given by coherent vectors 
multiplied by geometric (Berry) phase factors for a confined ion or for the center of mass 
(CM) of a system of identical ions. Thus we obtain an energy spectrum for the time 
independent Hamiltonians (the case of Penning type traps) or a quasienergy spectrum 
for time periodical Hamiltonians (in case of dynamical Paul traps). The coherent states 
are parametrized by the solutions of the classical equations of motion. 

These results have been extended to the case of an ion confined in a quadrupole 
electromagnetic trap with anharmonicities (a situation according to reality), in order 
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to study semiclassical dynamics associated to ions. Identification of equilibrium 
configurations is of interest for quantum information processing. We liave obtained 
tlie explicit expression of the quantum Hamiltonian associated to a charged particle, 
located in an nonlinear electromagnetic trap (with anharmonicities), as a function 
of the Lie algebra generators for the radial and axial symplectic groups. The 
Schrodinger equation solution was found for the Hamiltonian we considered. An 
analytical potential with cylindrical symmetry was chosen. The Hamiltonian was 
particularized to the case of combined, quadrupole and octopole traps. In case of the 
pseudopotential approximation, the minimum points of the dequantified Hamiltonian 
define the equilibrium configurations for trapped ions, which enable implementing and 
scaling of quantum logic for larger number of ions. 

Although the Schrodinger equation is nonlinear, the Hamiltonian (Husimi) function 
generally determines nonlinear equations of motion which can be investigated using 
the theory of differential dynamical systems. This theory enables qualitative global 
investigation of Hamilton families, depending on various control parameters. In 
particular it allows to describe classical chaos using Morse theory, bifucation theory 
and catastrophe theory. Thus quantum chaos is described using the methods of study 
for the deterministic classical chaos. This description is possible through implicit 
codification in the Husimi function of the spectral properties of the Schrodinger operator. 
Dequantization by means of coherent states enables to achieve exact connexions between 
symmetrical quantum and classical dynamical systems. For an ion confined within a 
quadrupole trap, the requantification corresponds to a geometrical quantification. 
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